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. black hole realized in nature. In the framework of acoustic geometry we study axially 

symmetric accretion and wind of a rotating astrophysical black hole or of a neutron 
!y— j ■ star assuming isentropic flow of a fluid described by a polytropic equation of state. In 

particular we analyze the causal structure of multitransonic configurations with two 
*q | sonic points and a shock. Retarded and advanced null curves clearly demonstrate the 

presence of the acoustic black hole at regular sonic points and of the white hole at the 
shock. We calculate the analogue surface gravity and the Hawking temperature for the 
inner and outer acoustic horizons. 
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Abstract 

Transonic accretion onto astrophysical objects is a unique example of analogue 



1 Introduction 



Analogue gravity has gained a lot of popularity in the last few years pQ (for a recent review 
and a comprehensive list of references, see (2j) for two obvious reasons: First, it provides a 
laboratory for testing some exotic effects of general relativity, such as black holes, wormholes, 
warped space-time, etc. Second, since the scale of analogue gravity is not fixed by the Planck 
length, there is hope that quantum gravity effects, which normally are important beyond 
the Planck scale only, may be seen on larger length scales. 

A particular example of analogue gravity is acoustic geometry in which an inhomogeneous 
fluid flow provides a curved space-time background for the propagation of sound waves. 
The most appealing aspect of acoustic geometry is a relatively easy way to construct an 
acoustic horizon which is an analogue of the black hole horizon in general relativity. In this 
way a number of interesting features which normally are extremely difficult to observe in 
astrophysical black holes, can be studied using an acoustic analogue of the black hole built 
in the laboratory. 

In all papers dealing with acoustic geometry the flows with desired properties are ar- 
tificially constructed. So far, these constructions have been only hypothetical since the 
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necessary conditions for the actual observation of the horizon and the measurements of the 
corresponding analogue Hawking temperature are still difficult to realize in the laboratory. 

Transonic accretion onto astrophysical black holes is a unique example of analogue grav- 
ity realized in nature [31 Ej. Spherical accretion was first studied in terms of relativistic 
acoustic metric by Moncrief [Sj. Recently, several aspects of analogue horizon radiation have 
been investigated for a spherically symmetric general relativistic |3] and post-Newtonian jE] 
accretion flow. In this paper we study acoustic geometry of an axially symmetric accretion 
disk round a rotating astrophysical black hole or round a neutron star. The flow considered 
in this paper is a relativistic generalization of the acoustic black hole analogues. To study 
acoustic geometry in most simple terms, we assume isentropic accretion of the fluid described 
by a polytropic equation of state. Compared with a spherically symmetric flow, an axially 
symmetric flow exhibits a more complex structure, which differs from the spherically sym- 
metric flow in two aspects. First, in axially symmetric accretion the acoustic horizon and the 
stationary limit surface generally do not coincide and second, multitransonic configurations 
with more than one acoustic horizons are possible. 

In section |21 we outline the basic relativistic hydrodynamics necessary for the study of 
accretion in terms of acoustic geometry. In section El we describe the relativistic equations 
governing the accretion flow. In section 0] we discuss the formation of a standing shock and 
acoustic horizons at sonic points and in section El we study the acoustic causal structure of 
multitransonic accretion and wind. Section |H1 is devoted to a calculation of the analogue 
surface gravity and in section HI we briefly discuss the acoustic superradiance. Finally, we 
conclude the paper with section |H1 

2 Preliminaries 

For the sake of completeness and consistency of the paper, as a starting point we review some 
well-known basic results in relativistic hydrodynamics and acoustics. We restrict attention 
to a stationary gravitating isentropic fluid in a stationary background metric g^ v . We work 
in units c = H = — &b — 1 and use the convention of positive-time-negative-space metric 
signature. 

2.1 Relativistic kinematics of the flow 

Let m m denote the 4- velocity field, i.e., the unit vector field tangent to the flow streamlines. 
In general, a 4- velocity can be expressed in terms of the 3- velocity components |Hj 







with the induced three-dimensional spatial metric 
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and the Lorentz factor 

^ 2 = I^> « 2 = 7««V; i,j = 1,2,3. (2.4) 
Let {Ej,} denote a set of hypersurfaces defined by 

w 2 _a 2 = 0, (2.5) 
where a is a constant, < a < c. Each is timelike since its normal 

rif, oc d^v 2 (2.6) 

is spacelike and may be normalized as n^n^ = —1. The 4- velocity at each point may be 
decomposed into normal and tangential components with respect to £„: 

u » = v± j ± n» + 7± L^, (2.7) 

where 

2/* = — ar + n"n>,; 7± = — L= . (2.8) 

t± v 1 ~ w i 

Next, we consider a two-dimensional axisymmetric flow in a stationary axisymmetric 
space-time. Physically, this may be a model for a black hole accretion disk, or for a draining 
flow in an axially symmetric bathtub. Using the notation v l = v r , v 2 = and t> 3 = v z = 0, 
the 4-velocity is given by 



u 



* = l(^-^v r ,vt,0), (2.9) 



<V = 7 ( 0rr«n -7= ~ (— ~ 94><t>) V ^ ) > ( 2 - 10 ) 

Vfi'oo #00 



w f = (^Mo - 9m) v I\ 7y = t - ^ 2 • ( 2 - 12 ) 



Then, the 3-velocity squared may be written as 

v 2 = v 2 - g rr v 2 , (2.11 

with 

?; 2 = (a 2 ^/ann — oj.j.)v 2 : ■y 2 = 

l-v 

For a hypersurface of constant v defined by ()2.5|) the decomposition ()2.7|) yields 

^l = -^7j; 7 2 = 7l7y- (2-13) 

It is useful to define a corotating frame as a coordinate frame of reference in which vn (or 
Vf/,) vanishes and, as a consequence v± = v. This frame is equivalent to the frame comoving 
with ZA In other words, in the corotating frame the vector L M takes a simple form 



The transition to a corotating frame is achieved by a simple coordinate transformation 

d(j) = d(j)' + Qdt , (2.15) 
where the angular velocity Q is defined as 

n = -. (2.16) 

More generally, we may define a coorbiting frame as a coordinate frame of reference 
comoving with ZA Hence an axially symmetric coorbiting frame is corotating. 



2.2 Acoustic geometry 

The acoustic metric tensor and its inverse are defined by |S] 



= — \9» v - (1 - c 2 s )u,u u ], GT = ^ 
wc* n 



- (1 - ^Kit" 



(2.17) 



where n is the particle number, w = {p + p)/n the specific enthalpy of the fluid and c s the 
speed of sound. In analogy to the event horizon in general relativity, the acoustic horizon 
may be defined as a time-like hypersurface whose normal is null with respect to the 
acoustic metric, i.e., 

G"%n„ = 0. (2.18) 
Equivalently, we may define the acoustic horizon as a hypersurface defined by the equation 

El 

v\ - c\ = 0. (2.19) 

The metric defined by (J2.17|) is very general and with no further assumptions the determi- 
nation of acoustic horizons may be nontrivial. However, the matter may greatly simplify if we 
assume a certain symmetry, such that the analogy with the familiar stationary axisymmetric 
black hole may be drawn. 

We consider the acoustic geometry that satisfies the following assumptions: 

1. The acoustic metric is stationary. 

2. The flow is symmetric under the transformation that generates a displacement on the 
hypersurface £„ of constant v along the projection of the flow velocity. 

3. The metric g^ is invariant under the above transformation. 

The last two statements are equivalent to saying that a displacement along the projection 
of the flow velocity on E v is an isometry. The axisymmetric flow considered in section 12.21 is 
an example that fulfils the above requirements. 

It may be easily shown that, for an axially symmetric flow, the metric discriminant 

^ = — GqqG^ (2.20) 

vanishes at the acoustic horizon. 
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In analogy to the Kerr black hole we may also define the acoustic ergo region as a region 
in space-time where the stationary Killing vector £ becomes space-like. The magnitude of 
defined with respect to the acoustic metric, is given by 

71 

lie ll= G^ec = —e^i 2 (c 2 s -v 2 ). (2.2i) 

wc s 

This becomes negative when the magnitude of the flow velocity v exceeds the speed of sound 
c s . The boundary of the ergo- region defined by the equation 

v 2 - c\ = 0. (2.22) 

is a hypersurface called stationary limit surface [9 j . Equivalently, the stationary limit surface 
may be defined as a hypersurface at which 

Goo = 0. (2.23) 

In the following sections we often use the terms subsonic, supersonic and transonic. These 
notions in relativistic hydrodynamics make sense only with respect to a specified frame of 
reference. For example, a flow which is supersonic with respect to an observer at rest is 
subsonic with respect to a comoving observer. In an axially symmetric flow it is convenient 
to use the corotating frame of reference, i.e., the frame in which vanishes, to define a locally 
subsonic or supersonic flow in the following way: the flow at a radial distance r is said to 
be subsonic or supersonic if v(r) < c s (r) or v(r) > c s (r), respectively. It may be shown 
that the above definition is consistent with another definition valid in any other reference 
frame that preserves temporal and axial isometry: the axisymmetric flow is said to be locally 
subsonic or supersonic if v±(r) < c s (r) or v±(r) > c s (r), respectively. Alternatively, as we 
demonstrate in sectional we may use equation ()2.20|) to define a flow with V > as subsonic 
and with T> < as supersonic. Finally, a flow is said to be transonic if a transition from a 
subsonic to a supersonic flow region, or vice versa, takes place. The point of transition is 
called the sonic point if the transition is continuous. A discontinuous supersonic to subsonic 
transition is called shock. 



2.3 Propagation of acoustic perturbations 

Propagation of acoustic waves in an inhomogeneous fluid moving in curved space-time is 
governed by the massless wave equation 



where 



^(V^GG^^O, (2.24) 



G = detG IMU . (2.25) 

It may be shown that equation (|2.24j) follows from the free wave equation that describes the 
propagation of acoustic perturbations in a homogeneous fluid in flat space-time 
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with a; = c s t, x l = (x,y,z) and r/^ = (1,-1,-1,-1). Equation (|2.24|) is obtained by 
replacing rj^ by G^ LU and d/dx^ by the covariant derivative V M associated with the metric 
Gfj, u . In this way we introduce the "acoustic equivalence principle" in full analogy with 
general relativity. 

We also define the stress tensor of the field if 

= d v <p - l -G afi d^ dpip G, lv (2.27) 

which satisfies G^V^T^ = 0. 

By standard arguments Elj it follows that the lines of acoustic propagation are null 
geodesies with respect to the acoustic metric. In other words, the wave vector which is 
timelike and may be normalized to unity, is null with respect to the acoustic metric, i.e., 

G^kW = 0, (2.28) 

and satisfies the acoustic geodesic equation 

FV M fc v = 0, (2.29) 



where 
and 



v„r = + r"k a (2.30) 



— T^G^^iG^a + G a p.^ — G^p). (2.31) 

However, as we have said, k 11 is timelike with respect to the space-time metric, so acoustic 
perturbations propagate along timelike curves. 

Next we show that the perturbations that propagate parallel to the streamlines cross 
the horizon and those that propagate in the opposite direction are trapped. We assume 
that the space-time metric is stationary and the flow is steady. At any fixed time to the 
flow congruence crosses the two-dimensional surfaces {t> 2 =const, t = to}. In particular, the 
acoustic horizon defined by v\ = c 2 at fixed t = t is one of these surfaces. The hypersurface 
v\ = c 2 s is of course timelike since its normal d^v 2 ^ — c 2 ) is a spacelike vector. Let 
be a timelike unit vector tangent to the lines of acoustic propagation. Then k^ may be 
decomposed in a similar way as u M in equation (|2.7|) . 

p = klk7l n + 7fcPj (2.32) 

where 

P=V + ^%; Tfc (2.33) 
Ik VI - k 2 

Now we restrict attention to the perturbations for which = ZA In the coorbiting frame 

defined in section 12 .2\ these perturbations propagate either parallel or antiparallel to the 

streamlines, i.e., the spatial part of k^ is parallel or antiparallel to the 3- velocity of the 

fluid. As a consequence of ()2.28|) we have k^u^ = 7 Cs . This equation states the fact that 

acoustic perturbations propagate with the velocity of sound c s with respect to the fluid. At 

the horizon, as may be easily shown, k = 2c s /(l + c 2 ) for the perturbations going parallel, 

and k = for those going antiparallel to the streamlines. Since in the first case we have 

k^n^ ^ 0, the acoustic perturbations that propagate parallel to the streamlines will cross 

the horizon. The perturbations going antiparallel will satisfy k^n^ = at the horizon and 

hence will be trapped. 
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2.4 Quantization of phonons and the Hawking effect 

The purpose of this section is to demonstrate how the quantization of phonons in the pres- 
ence of the acoustic horizon yields acoustic Hawking radiation. The acoustic perturbations 
considered here are classical sound waves or phonons that satisfy the massles wave equa- 
tion (|2.24jl in curved background with the metric given by (j2.17|) . Irrespective of the 
underlying microscopic structure acoustic perturbations are quantized. A precise quantiza- 
tion scheme for an analogue gravity system may be rather involved jTTj. However, at the 
scales larger than the atomic scales below which a perfect fluid description breaks down, the 
atomic substructure may be neglected and the field may be considered elementary. Hence, 
the quantization proceeds in the same way as in the case of a scalar field in curved space 
[T2*] with a suitable UV cutoff for the scales below a typical atomic size of a few A. 

For our purpose, the most convenient quantization prescription is the Euclidean path 
integral formulation. Consider a 2+1-dimensional disc geometry. The equation of motion 
(J2.24)) follows from the variational principle applied to the action functional 



S[<p] = / dtdrd^V^GG^d^dufp. 



(2.34) 



We define the functional integral 



Vipe 



(2.35) 



where Se is the Euclidean action obtained from (|2.34j) by setting t — ir and continuing the 
Euclidean time r from imaginary to real values. For a field theory at zero temperature, the 
integral over r extends up to infinity. Here, owing to the presence of the acoustic horizon, 
the integral over r will be cut at the inverse Hawking temperature 2ir/n where k denotes 
the analogue surface gravity. To illustrate how this happens, consider, for simplicity, a 
nonrotating fluid (v<f, = 0) in an arbitrary static, spherically symmetric space-time. It may 
be easily shown that the acoustic metric takes the form 



ds 2 



l-v 



y2 1-C 2 

g 00 dt 2 - 2v- W-9oo9rr drdt + 



c 2 v 2 



l-v 2 



l-v 2 



-9r 



dr 2 



(2.36) 



where the background metric coefficients g w are functions of r only, and we have omited the 
irrelevant conformal factor n/(wc s ). Using the coordinate transformation 



dt — > dt + 



v(l-c 2 ) 



9rr , 

dr 

9oo 



(2.37) 



we remove the off-diagonal part from ()2.36|) and obtain 



ds 2 



l-v 2 



9oo dt 2 + 



-2\2 



2 2 , 
C S V + 



9r 



dr — r' L 



Next, we evaluate the metric near the acoustic horizon at r 
r — rh at first order 

d_ 

dr 



ci-v z w 2c s — (c. 



(2.38) 

rh using the expansion in 
r - r h ) (2.39) 
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and making the substitution 



l id, _ , 
r Th ~ 2c s (l-c*)g rr dr [Cs V > 



p\ (2.40) 



where p denotes a new radial variable. Neglecting the subdominant terms in the square 
brackets in ()2.38|) and setting t = ir, we obtain the Euclidean metric in the form 

ds% = n 2 p 2 dr 2 + dp 2 + rldcf 2 , (2.41) 

where 

K= T ^ x f^Uc.-v) . (2.42) 
1 - eg V 9rr 9r h 

Hence, the metric near r = is the product of the metric on S 1 and the Euclidean Rindler 
space-time 

dsl = dp 2 + p 2 d(Kr) 2 . (2.43) 

With the periodic identification r = r + 2-k/k, the metric (J2.43JI describes M 2 in plane polar 
coordinates. 

Furthermore, making the substitutions p = e KX /n and <fi = y/r^ + Tc, the Euclidean action 
takes the form of the 2+1- dimensional free scalar field action at nonzero temperature 

/•27t/k poo poo 1 

SeM = / dr dx dy-(d^) 2 , (2.44) 

JO J-oo J-oo ^ 

where we have set the upper and lower bounds of the integral over dy to +oo and — oo, 
respectively, assuming that is sufficiently large. Hence, the functional integral Z in (j2.35j) 
is evaluated over the fields ip(x,y,r) that are periodic in r with period 2tt/k. In this way, 
the functional Z is just the partition function for a grandcanonical ensemble of free massless 
bosons at the Hawking temperature Th = k/(2tt). However, the radiation spectrum will not 
be exactly thermal since we have to cut off the scales below the atomic scale ^H] • The choice 
of the cutoff and the deviation of the acoustic radiation spectrum from the thermal spectrum 
is closely related to the so-called transplanckian problem of Hawking radiation [T4] . 

2.5 Conservation laws 

Our basic equations are the continuity equation 

(mz% = -L^(v^n/) = (2.45) 



-9 

and the energy-momentum conservation 

T» v . v = 0. (2.46) 
Applied to a perfect, isentropic fluid, this equation yields the relativistic Euler equation [T3] 

u u (w Ufl ). >v - d^w = 0. (2.47) 
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The quantities c s , and may be discontinuous at a hypersurface E. In this case, 
equations (|2.45jl and ()2.46|) imply 

[frm"]] n„ = 0, (2.48) 

[[Tn]^ = 0, (2.49) 
where [[/]] denotes the discontinuity of / across E, i.e., 

[[/]]=/2-/i, (2-50) 

with / 2 and /i being the boundary values on the two sides of E. For a perfect fluid, equations 
(|2~4*%)) and (|2~4^)) may be written as 

[KltJ] = 0, (2.51) 

[[ToXW = [[(P + P)«o«±7±]] = 0, (2.52) 
[P^n"]] = [[(p + pKt! +p]] = 0. (2.53) 
Equations ()2.51J) - (j2.53[) are referred to as the relativistic Rankine-Hugoniot equations [16J. 



2.6 Constants of motion 

The acoustic metric and the background space-time metric generally do not share the isome- 
tries. For example, the axial symmetry of a tilted accretion disc round a rotating black hole 
obviously does not coincide with that of the black hole. 

Here we assume that the flow and the space time metric share a specific symmetry 
generated by a Killing vector x M . 

Proposition 1 Let % M be a Killing vector field and let the flow be isentropic and invariant 
under the group of transformations generated by x^ ■ Then the quantity wx^Ufi is constant 
along the stream line. 

Proof: Multiplying the Euler equation (|2.47j) by % M we obtain 

u u d v {wx»Ur) - uWwXw - X^w = 0. (2.54) 

The second term in ()2.54j) vanishes by Killing's equation and the last term vanishes by 
symmetry. Hence 

u u d v (w X ! i u ll )=0, (2.55) 

as was to be proved. □ 

Remark Since, by assumption, the flow shares the symmetry of g^ u , the Killing vector 
X M is also Killing with respect to the acoustic metric. 
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2.7 Polytropic gas 

Consider a gas of particles of mass m described by the polytropic equation of state 

p = m/Cn r , (2.56) 

where fC is a constant which, as we shall shortly see, is related to the specific entropy of the 
gas. Using this equation and the thermodynamic identity 

dw = Td (-) + -dp, (2.57) 
\n/ n 

the energy density p and the entropy density s may also be expressed in terms of n. For an 
isentropic flow, it follows 

p = mn ^1 + f^-J nF ~^j ■ ( 2 - 58 ) 

For an adiabatic process, the entropy density is proportional to the particle number 
density, i.e., s = an. Without further assumptions, the entropy density per particle a is an 
undetermined constant. If, in addition to ()2.56|) . the Clapeyron equation for an ideal gas 
holds 

p = nT, (2.59) 
and assuming that 7 is temperature independent, then the entropy per particle is given by 

a = ^ log(m/C) + C + ^ , (2.60) 

where ( is a constant that depends only on the chemical composition of the gas. Hence, the 
quantity /C measures the specific entropy. 

It is convenient to express all relevant thermodynamic quantities in terms of the adiabatic 
speed of sound defined by 

2 n dw 



s w dn 



(2.61) 



We find 

w r — 1 . 

m 1 — 1 — Cg 

1 /T-l V 1 ^ 
n = r)Krr=i _ ij , (2.63) 

1 /r — 1 \ 

p = rmfK,-—^ — 1 , (2.64) 



where 



V = ( ^A) r " • (2-65) 
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3 Accretion disk 



For a flow of matter with a non-zero angular momentum density, spherical symmetry is 
broken and accretion phenomena are studied employing an axisymmetric configuration. Ac- 
creting matter is thrown into circular orbits, leading to the formation of accretion disks 
around astrophysical black holes or neutron stars. The pioneering contribution to study 
the properties of general relativistic accretion disks is attributed to two classic papers by 
Bardeen, Press and Teukolsky [TH] and Novikov and Thorne 

If the accreting material is assumed to be at rest far from the accreting black hole, the flow 
must exhibit transonic behaviour since the velocity of fluid particles necessarily approaches 
the speed of light as they approach the black-hole event horizon. For certain values of the 
intrinsic angular momentum density of the accreting material, the number of sonic points, 
unlike in spherical accretion, may exceed one, and accretion is called multi-transonic. Study 
of such multi-transonic flows was initiated by Abramowicz and Zurek [20J. Subsequently, 
multi-transonic phenomena in black hole accretion disks and formation of shock waves as 
a consequence of these phenomena, have been studied by a number of authors (for further 
details and for references see, e.g., [H l2"Tl l2"2l 123J). 



3.1 Axisymmetric disc accretion 

In this section, we outline the formalism for calculating the location of sonic points as a 
function of the specific flow energy S, the specific angular momentum A and the polytropic 
index T [221 E3] applied to axisymmetric accretion. We assume that the accretion flow is 
described by a disk placed in the equatorial plane of a rotating astrophysical object, e.g., 
a Kerr black hole or a neutron star, of mass M. Specifying the metric to be stationary 
and axially symmetric, the two generators £ M and M of the temporal and axial isometry, 
respectively, satisfy equation (|2.55|) and hence the quantities £ = wuq and L = wu^ are 
constant along the streamlines. Given the equation of state of the fluid, the constants of 
motion 8 and L fully determine the accretion flow. It is convenient to introduce another 
constant 

A = = (3.1) 

5 Uq 

which is called specific angular momentum. It may be easily shown that the angular velocity 
fl defined by f!2.16|) is expressed in terms of A as 

Q = _go0 + Agoo (32) 



With the help of this equation, either from the normalization condition u^u^ = 1 or directly 
from (l2~ToT) with (j2~T2>(j2~T3|) . we find 

«o = 7i|, (3-3) 
where, using the notation of JH1 ED] , A and B expressed in Boyer-Lindquist coordinates are 

A = 9o4, - 9009M = r 2 - 2Mr + a 2 , (3.4) 
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r, / , N / 2 2 2Ma 2 \ Ma / 2M\ . . 

5 = -{Qh + 2A^ + A 2 (?oo) = ( r 2 + a 2 + —j - 4A— - A f 1 - —J . (3.5) 

It is useful to introduce the particle accretion rate N which, owing to the particle number 
conservation, does not depend on the radial coordinate r. Besides assuming axial symmetry, 
we also assume that the disk thickness h is constant and the distribution of matter in the disk 
does not depend on the ^-coordinate along the rotation axis 1 . We start from the expression 

N = J nu^dE^ (3.6) 

for the number of particles the worldlines of which cross a hypersurface S. Defining the 
hypersurface at rest as a cylinder of radius r and height h and by making use of the axial 
symmetry and the continuity equation ()2.45j) integrated over the disk volume, we find 

N/h = -2it^g~nu r . (3.7) 

Using ()2.13|) we obtain 

N/h = -27rVAnu ± 7i. . (3.8) 

Thus, the two conditions 

N/h = -2irnv ± ^±A 1/2 = const , (3.9) 
A 1 / 2 

£ = wuq = u>7_l —, , = const , (3.10) 

together with the equation of state, fully determine the kinematics of the accretion disk. In 
other words, given the equation of state, e.g., in the form of (|2.56|h the metric and a set of 
parameters {A, £, N/h}, the flow variables v± and c s may be determined as functions of r. 



3.2 Non-axisymmetric accretion 

In this work, we concentrate only on axisymmetric fluid accretion, for which the orbital 
angular momentum of the entire disc plane remains aligned with the angular momentum 
of the compact object at the origin. In a strongly coupled binary system (with a compact 
object as one of the components), accretion may experience a non-axisymmetric potential 
because the secondary donor star exerts a tidal force on the accretion disc round the compact 
primary. In general, a non-axisymmetric tilted disc may form if the accretion takes place 
out of the symmetry plane of the spinning compact object. The matter in such a misaligned 
disc will experience a torque due to the general relativistic Lense-Thirring effect |24j . leading 
to the precession of the inner disc plane. The differential precession may cause stress and 
dissipative effects in the disc. If the torque remains strong enough compared with the internal 
viscous force, the inner region of the initially tilted disc may be forced to realign itself with 
the symmetry plane of the central accretor. This phenomenon of partial realignment, out to a 
certain radial distance known as the transition radius or the alignment radius, of the initially 

1 A physically more realistic calculation would involve solving the relativistic Euler equation (|2.47() in the 
vertical direction [23]. This work is in progress and will be presented elsewhere. 
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non-axisymmetric disc is known as the "Bardeen-Petterson effect" The transition radius 
can be obtained by balancing the precession and the inward drift or the viscous time scale. 

An astrophysical accretion disc subject to the Bardeen-Petterson effect becomes twisted 
or warped. A large scale twist in the disc modifies the emergent spectrum and can influence 
the direction of the quasar and micro-quasar jets emanating out from the inner region of the 
accretion disc 

A twisted disc may be thought of as an ensemble of rings of increasing radii, for which the 
variation of the direction of the orbital angular momentum occurs smoothly while crossing the 
alignment radius. A system of equations describing such a twisted disc has been formulated 
\27\ . and the time scale required for a Kerr black hole to align its angular momentum with 
that of the initially misaligned accretion disc, has been estimated f2E]|. Numerical simulations 
using the three-dimensional Newtonian smoothed particle hydrodynamics (SPH) code [29, 
as well as a fully relativistic framework [3U| reveal the geometric structure of twisted discs. 

Note, however, that as long as the acoustic horizon forms at a radial length scale smaller 
than that of the alignment radius, typically of order 100 - 1000 Schwarzschild radii according 
to the original estimate of Bardeen and Petterson [22], one need not implement the non- 
axisymmetric geometry to study the analogue effects. 

4 Acoustic horizon 

The flow considered here may be regarded as a relativistic generalization of the two-dimen- 
sional draining bathtub flow with a sink at the origin. Instead of a sink, at the origin we 
have an astrophysical black hole or a neutron star playing the role of the sink. Hence, our 
accretion flow is a general-relativistic black-hole analogue. 

In the axially symmetric flow, the acoustic horizon is a stationary circle of radius r^ 
where, according to (|2.19p . the radial component of the flow experiences a transition from a 
subsonic to a supersonic region. As we shall shortly demonstrate, in our perfect polytropic 
gas model such a transition may be realized in three ways: 

Smooth transition Both c s and v± and their derivatives are regular at the transition point. 
In this case, the transition point is called regular sonic point or simply sonic point. The 
points known in the literature as saddle-type or X-type critical points belong to this 
class. The isolated sonic points known as centre-type or O-type also belong to this 
class. 

Singular transition Both c s and v± are continuous but their derivatives diverge. The 
transition point is called singular sonic point. 

Discontinuity Both c s and v± are discontinuous as a consequence of discontinuity in the 
equation of state. The point of discontinuity is called shock. 

The second and the third case are consequences of the idealization of assuming an irro- 
tational inviscid perfect fluid. In reality, the divergences will be absent because one or more 
simplifications become invalid as v± approaches c s or at a shock. However, the mathematical 
divergences indicate that the corresponding divergent physical quantities, such as analogue 
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5 10 50 100 500 1000 r 

Figure 1: Mach number v±/c s versus r (in units of the Schwarzschild radius) for various 
N/h. The entropy measure /C corresponding to the flows marked by (a), (b) and (c) equals 
0.0181289, 0.0188906. 0.0213815, respectively. 




surface gravity, may in a real physical situation be much larger than previously expected 

To find solutions with at least one regular sonic point, we take the derivative of (J3.9)) and 
(|3.10|) with respect to r, which yields 

/ . \ / 

AT I AT \ 

= 0, (4.1) 

w'u + w(u )' = 0. (4.2) 

Substituting N/n and Uo from (J3.9J) and (|3.10|) and using ()2.6H) . equations ()4.1|) and ()4.2|) 

may be written as a system of linear algebraic equations for w'/w and v'_ L /v± 

- + chl— = -c 2 s --r , (4.3) 
w v± 2 A 

hw,7| — = — H . (4.4) 

w ± ± v ± 2 A 2B K ' 



Obviously, the determinant of the system vanishes at the acoustic horizon owing to ()2.19|) . 
The system will have a regular solution at the horizon if and only if the right-hand sides of 
the two equations become equal at the horizon. In this case, the algebraic system ()4.3|) - (j4.4j) 
degenerates at the horizon into one equation. Hence, the condition 

(l-c^-f = (4.5) 
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must hold at the horizon. If this condition is not satisfied, any solution to (|4.3|) - (|4.4|) will 
have no acoustic horizon at all or will have irregular acoustic horizons, i.e., such that the 
derivatives v' ± and c' s approach infinity as one approaches the horizon. 

The regularity condition ()4.5|) gives c s at the horizon in terms of r^ but the actual value 
of ?"h is yet unknown. To find r^, we need to specify the equation of state. Given the equation 
of state, we can express w in terms of c s , substitute c s from ()4.5|) into (j3.10J) and thus obtain 
an equation for r. For example, equation ()2.62)1 . which follows from the polytropic equation 
of state (|2.5fjjl . yields 

A solution to this equation for fixed A and £ determines the horizon radius r^. Equation (|4.fij) 
may, in general, have more than one solution. However, not all real solutions are physically 
acceptable. The positivity of in (j3.3|) requires 

g H + 2\g Q4> + A 2 #oo < 0. (4.7) 

For a given A, this inequality fixes the range of acceptable r. 

Once we have and c s at the horizon, the constant N/h may also be fixed from equa- 
tion ()3.9|) evaluated at the horizon. Equations ()3.9|) and (|3.10|) then determine the radial 
dependence of c s and v±. It is important to emphasize that this procedure determines only a 
subset of solutions to equations ()3.9|) and ()3.1U|) (or equivalently ()4.3|) and (|4.4|) ) which pos- 
sess sonic points that satisfy the condition (|4.5|) . Any other solution will have singular sonic 
points or no sonic point at all. This may be seen in figures Hand El where two typical sets of 
solutions with different topology depending critically on A are plotted for a fixed £. Solutions 
are expressed in terms of the so-called Mach number defined as the ratio M. = vj_/c s . The 
curves crossing the M. — 1 line are transonic and the crossing points are the location of the 



15 



M 




t 1 1 T — I — i — r~| — 

r r 10 

hi shl 



- 1 1 1 1 1 — t — I — r-| 1 1 1 1 r- 

r r 100 r 

sh2 h2 -l 



Figure 3: Transonic accretion with two sonic points at r^i and r\a and a shock in between 
at ^sh2 (solid line). The dashed line represents an alternative flow with a (unstable) shock 
at r s hi . The parameters are as in figure [T] 



acoustic horizon. The plots also include the configurations that do not satisfy ()4.5)l at the 
acoustic horizon. The plots describe both wind and accretion, depending on the sign of N 
or equivalently of v±. Each line on the plot represents wind if N < (v± > 0) or accretion if 
N > (vj_ < 0). The flow described by line (a) has an X-type sonic point. The flow marked 
by (b) has two branches, each passing an X-type and ending at a singular sonic point. An 
isolated O-type sonic point is marked by (c). 

Obviously, a smooth flow will have at most one regular sonic point.. Hence, a smooth 
multitransonic flow is not possible. However, a multitransonic flow with two sonic points 
and a shock in between may be realized if discontinuous transitions between different lines 
are allowed. 

The lines in the plots in figures ^ and El in general correspond to different constants N/h 
and if the entropy is conserved, there is no way to jump from one line to another. However, 
instead of keeping the constant K fixed, we fix N/h and allow /C to vary from line to line. In 
this way the jumps between the lines labelled by different JC will be allowed at those points 
where the Rankine-Hugoniot conditions are satisfied. 

The flows with singular sonic points, i.e., the points where v' ± and c' s diverge, will formally 
have acoustic horizons with infinite surface gravity. This property is generic for inviscid flows 
|3T] and the problem will not appear in a realistic fluid with viscosity. In our case, owing 
to the shock formation, the flow will actually never reach these points. We assume that a 
stationary accretion flow goes from infinity where it is subsonic to the first unstable orbit of 
the black hole. In the example depicted in figure ^ the flow represented by line (a) passes 
through a regular sonic point. If the Rankine-Hugoniot conditions are met, a jump from a 
supersonic to a subsonic regime will take place, provided the entropy of the subsonic flow is 
larger. Since K.u\ < /C( b ) < JCm, the entropy per particle a, which increases with K according 
to ()2.60|) . increases going from line (a) towards the point denoted by (c) and hence, a jump 
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Figure 4: Transonic wind with two sonic points at r^i and and a shock in between at 
f B hi (solid line). The dashed line represents an alternative flow with a shock at r s h2- The 
parameters are as in figure El 



from (a) to (b) or to any of the closed loops inside loop (b) is possible. However, the only 
transition that will eventually result in a stable stationary flow will be a jump from (a) to 
the subsonic branch of (b). A flow with a similar transition from (a) to one of the closed 
loops will inevitably end up at a singular sonic point, and hence cannot remain stationary. 

The exact location of the shock r s h for a particular flow configuration may be calcu- 
lated using the Rankine-Hugoniot conditions (|2.51|) - ([2.53)) assuming that the outer flow is 
described by (a) and the inner flow by the subsonic branch of (b). The first two conditions 
(I2.51|) and (|2.52j) are trivially satisfied owing to equations (13 .9|) and (|3.1U|) . The particle 
number conservation ()3.9|) across the shock may be written as 

JC'^Afx = K~ ^W 2 , (4.8) 

where the subscripts 1 and 2 denote, respectively, the inner and outer boundary values at 
the shock and 

M = V ( — " 1 ) v ±1± (4.9) 

is a kinematic quantity the change of which across the shock depends on the change of v± and 
c s only. Hence, to make the particle number conserved, a discontinuity of N is compensated 
by a corresponding discontinuity of K such that equation ()4.8|) holds. As a consequence, the 
entropy per particle a being a continuous function of K, is not conserved across the shock. 
The third Rankine-Hugoniot condition ()2.53|) may now be written as 
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v r+1 ( ^ - if" 1 + (r- 1)4 ( ^ - i) r_1 ^ 



2 J 1 V "I"** \ r 2 



K.2 



(4.10) 



and may be solved for r. In figures El and 0] we show the solutions for the accretion and 
wind, respectively, corresponding to the flows depicted in figures [T] and El The two possible 
locations of the shock are shown by vertical lines at r s hi and r s ^ 2 with arrows indicating the 
flow orientation. According to a standard local stability analysis jS21; one can show for a 
multitransonic accretion that only the shock formed in between the middle (indicated by (c) 
in figures and EJ) and the outer sonic point is stable. Hence, in our accretion example in 
figure El the shock at r s h2 is stable and the one at r s hi is unstable. 

5 Acoustic metric and causal structure 

The geometric aspects of the axially symmetric accretion are conveniently described in terms 
of the acoustic metric defined by ()2.17|) and the metric discriminant T> defined by ([2.20)1 . 
In figure El we plot the components of the metric tensor as functions of r calculated using 
equation (|2.17J) for the transonic accretion flow represented by the solid line in figure El and 
similarly in El for the transonic wind of figure |U For the same flow configurations depicted 
in figures El and |U we also plot V as a function of r in figures [7| and El respectively. Clearly, 
T> is positive (negative) in the subsonic (supersonic) region. Acoustic horizons are located 
at the points where the discriminant changes its sign. Depending on the slope of T> at 
these points passing in the direction of the flow, we call them acoustic black hole (if the 
slope is descending) or white hole (if the slope is ascending) horizons. Indeed, for both the 
multitransonic accretion (figure EJ) and the wind (figure EJ), the flow at the inner and the 
outer sonic points denoted by r hl and r h2 , respectively, goes from subsonic to supersonic; 
so, according to the classification of Barcelo et al [33], there we have acoustic black hole 
horizons. At the discontinuity r shl or r sh2 the flow goes from supersonic to subsonic and 
hence observers in the subsonic region of the accretion (between r^i and r s h2 in figures El and 
0), and of the wind (between r shl and r h2 in figures E| and EJ) experience a white hole. 

In figures [7| and El we also plot the metric component Goo as a function of r (represented 
by the dashed line). The point r where Goo changes its sign is the radius of the stationary 
limit surface and the boundary of the acoustic ergo region. 

In order to study the causal structure of the axially symmetric acoustic geometry, let 
us introduce retarded and advanced radial null coordinates u and w, respectively. We first 
make a transformation to the acoustic corotating frame by eliminating the term G^dtdcf) 
from the acoustic metric. This is achieved by transforming to 0' as 

d^' = d^ + ^dt, (5.1) 
G^ 

with the remaining coordinates unchanged. This yields the acoustic metric 

ds 2 = G' flu dx^dx u , (5.2) 

where 



Gqo = — 7; — ; G 0lj) = 0, (5.3) 
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Figure 7: Acoustic metric functions V (solid line) and G o (dashed line) versus r for the 
multitransonic accretion as in figure El with a shock at r s h2- Acoustic horizons and their type 
(from the point of view of an observer in the subsonic region) are indicated by BH (black 
hole) and WH (white hole). 




Figure 8: Acoustic metric functions T> (solid line) and Goo (dashed line) versus r for the 
multitransonic wind as in figure 0] with a shock at r s hi. 
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Figure 9: Acoustic causal structure of the transonic accretion of figureElwith the null curves 
u = const (solid lines) and w = const (dashed lines). 
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G Qr = G r0 = Go r — G r $— — , (5.4) 

with all remaining components of the acoustic metric unchanged. We then make the trans- 
formation 

du — dtA dr , (5.5) 

c+ 

dw — dt H dr , (5.6) 

where 



_ *fOr ^ V IfOr ^PO^rr /_ -s 

C± — ^ . {O.I ) 



The acoustic metric at constant 0' is then given by 

V 

ds 2 \,, — — — — dudw . (5.8) 

I <p'=const v 7 

The functions £ = £(r) at constant u or constant w, obtained by integrating 

t = - —dr, (5.9) 
J c± 



represent, respectively, the right-moving and left-moving sound rays or null-rays 

In figures El and ^1 we plot the u = const and w = const lines for the accretion of figure 
EJand the wind of figure 01 The divergence of lines at the acoustic horizons towards t = — oo 
demonstrate an acoustic black hole. Since the transition from a supersonic to a subsonic 
region is discontinuous at the shock, the lines have a cusp pointing in the direction of positive 
t. If the transition at the shock were smooth, the null lines would diverge towards t = +oo 
as in the white-hole example of Barcelo et al [HH]- Hence, the cusp indicates a white hole. 



6 Surface gravity and Hawking temperature 

One of the most interesting aspects of acoustic horizons is the existence of a surface gravity 
and the associated Hawking radiation. To calculate the analogue surface gravity, we need 
the derivatives c' s and v' ± at the acoustic horizon. Let us consider a regular sonic point, e.g., 
at rhi or rh2 in figures OB first. The case of shock will be discussed afterwards. To calculate 
c' s and v' ± at a regular sonic point, we have to solve two equations. The first one is obtained 
combining ()4.3|) and the derivative of ()2.62|) : 

4ci 2v' ± A' 

+ --, o=-°s^r- (6-1) 



r-i- c 2 i-c* a 

The second equation is obtained by eliminating w'/w from ()4.3j) and ()4.4j) and taking the 
derivative with respect to r of the equation thus obtained at the horizon. We find 

^l(c' s - v' ± ) + 2c s c' s ^ = (1 - cl) - ^ j - §■ + . (6.2) 
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Figure 11: Analogue Hawking temperature at the inner (left) and outer (right) acoustic 
horizon versus the specific angular momentum of the black hole. 



Equations (j6.1|) and ()6.2j) can now be easily solved for c' s and v' ± . 

The analogue surface gravity k may be calculated with help of the Killing field x M that 
is null on the horizon. We start from the expression [HI ED] 

G^n v ^-{G apX a X P ) = 2«X M , (6.3) 

where d/dn = n^d^ denotes the normal derivative at the horizon. Since the definition of 
surface gravity by equation (|6.3|) is conformally invariant in the calculations that follow 
we drop the conformal factor n/(wc s ) in G^ v . 

Consider the vector field defined by (|2.8j) . Its magnitude in the acoustic metric is 
given by 

\\L\\>=GrWL" = [l-{l-<%^=rttf-v* ± ). (6.4) 

Hence, the vector is null on the horizon. We now construct a Killing vector x^ m the 
form 

x? = e + n«if, (6.5) 

where the Killing vectors £ M and M are the generators of the temporal and axial isometry 
group, respectively. It may be easily verified that if the constant coincides with the 
angular velocity (|2.16|) evaluated at the acoustic horizon, then x M becomes parallel to L M 
and hence null on the horizon. Using x^ = y /x u XuL fJ- and the definition ()2.8j) we find 

1 1 y M 
G^n v = n" = £== (6.6) 

at the horizon. Equation ()6.3|h together with ()6.4j) and (|6.6|) yields 



" = r^V^ (< - 0> (6 ' 7) 
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Figure 12: Analogue Hawking temperature at the inner (left) and outer (right) acoustic 
horizon versus the specific angular momentum of the disk for the two topologies depicted in 
figure H] (solid line) and figure 121 (dashed line). 



where it is understood that the derivatives and the quantity x^X^ are to be taken at the 
horizon. The corresponding Hawking temperature Tah = l K l/(2vr) represents the temper- 
ature as measured by an observer near infinity. For A = a = 0, equation (jfi.7|) reduces to 
1)2.42)1 derived in a different way for a nonrotating fluid in a static background metric. 
The norm of the Killing vector x^ is easily calculated from (J6.5)) using ()3.2|) . We find 

V X ^ ~ r 3 + a 2 r + 2Ma? - 2\Ma 

We plot the ratio of the analogue Hawking temperature Tah to the black hole Hawking 
temperature Th at the inner and outer horizons as a function of the black-hole specific angular 
momentum a (figure ITT]) and of the disk specific angular momentum A ( figure IT2*j) . The solid 
and dashed lines represent the two topologies of figures UJ and EJ respectively. However, 
whereas the inner and the outer acoustic horizons for the topology of figure belong to the 
same multitransonic flow connected by a shock (the typical configuration depicted in figure 
EJ), the horizons for the topology of figure El belong to disconnected transonic flows marked by 
(a) and (b) in figure El Amazingly, the analogue surface gravity is smooth at the topological 
transition point. 

Note that the analogue Hawking temperature Tah of the outer horizon is much lower than 
the temperature Tah of the inner horizon, which in turn is, although close in magnitude, 
typically lower than the black hole Hawking temperature Th- This is due to the relatively 
low specific energy £ = 1.01, barely exceeding the rest mass. For a range of larger values of 
£ , e.g., close to 2, the temperature Tah would exceed Th as in the spherically symmetric 
accretion [3J Ej. However, we have found that in this range of large specific energies a 
multitransonic behaviour is absent. 

Calculation of the surface gravity at the shock is not possible in the present model with 
no viscosity. The shock may be viewed as a boundary separating two phases of the fluid. 
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The equation of state, and hence the sound speed, as well as the radial component of the 
fluid velocity, in each phase are different and exhibit discontinuity at the boundary. As a 
consequence, the surface gravity and the associated Hawking temperature at the shock will 
be formally infinite. However, in a more realistic situation, the phase boundary, owing to 
viscosity, will have a finite thickness and discontinuities will be smoothed out. The gradients 
of v± and c s will actually be finite but possibly very large. A similar situation occurs in a 
non-relativistic acoustic setup in which the boundary of a stable-phase bubble propagates in 
the surrounding metastable phase |3S] and in the 3 He-A superfluid where a superluminally 
moving domain wall soliton provides an analogue event horizon [HE1 EZ| ■ 

7 Acoustic superradiance 

The existence of an ergo region in the transonic accretion flow as described in sections 0] and 
El suggests the possibility of acoustic superradiance. Acoustic waves entering the ergo region 
and reflecting from the outer acoustic horizon may experience a coefficient of reflection 
greater than 1. This effect is a close analogue of superradiant scattering [SHI EH] which 
allows energy extraction from Kerr black holes. Acoustic superradiance has recently been 
extensively discussed in the context of the draining bathtub [30j Ej and its application to 
the Bose-Einstein condensate [32] . 

The quantitative investigation of acoustic superradiance in the accretion disc geometry 
may be performed in the usual way: First, the variables in (J2.24J1 may be separated using a 
cylindrical wave ansatz 

<p(t,r,(f>) = $(r)e-^V m *, (7.1) 

where u> is real and positive and m a positive integer. Equation (|2.24j) then reduces to 
an ordinary second-order differential equation for $(r) which may be solved numerically 
for suitable boundary conditions. However, this calculation may be quite involved for the 
transonic flow described in sections |U and El and would go beyond the scope of this paper. 
Instead, we would like to demonstrate here that the acoustic superradiance occurs in a 
relativistic transonic axisymmetric flow if 

< uj < mVL h , (7.2) 

where flh is the flow angular velocity (|2.16|) evaluated at the acoustic horizon. To show this, 
we follow the procedure outlined in [TO] . 

Using the stress tensor (|2.27|) and the stationary Killing vector £ M we first construct the 
analogue "energy current" 

Jv = T fa ,e", J» = GTJ V , (7.3) 

which is conserved, i.e. which satisfies V M J M = 0, V M being the covariant derivative associ- 
ated with G^ v . Hence if we integrate V M J M over the region S with the boundary dS, we find 
by Gauss's law 

/ dS^= [ dS^+ [ dS ll J"+ [ dS^+ [ dS^ = Q, (7.4) 

JdS JC JEi JS 2 J AH 
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where we have chosen dS such that it consists of the cylinder C of large radius R — > oo and 
height At, the disc Si at the bottom of the cylinder extended from the outer acoustic horizon 
radius to R, the time translated disk S 2 at the top, and the part ATC of the outer acoustic 
horizon 7i between Si and S 2 . The integrals over Si and S 2 cancel by time symmetry, the 
integral over C represents the net energy flow out of S to infinity, i.e., the difference between 
the outgoing and incoming energies E ont — E- in during the time At, whereas the integral over 
ATi, represents the net energy flow into the acoustic black hole. Thus we have 

£out = dS^ = [ dAdt X % , (7.5) 

J AH J 

where we have used = —G^ v x u as the outward directed normal to the horizon and x M is 
given by (jfi.5)l . The integral in ()7.5j) is taken over the proper volume element dAdt on the 
acoustic horizon. Using G^x^^ = which follows directly from (jfi.fij) . a straightforward 
algebra yields 

x% = Cd v <p = (dt<p + n h d^)dt<p. (7.6) 

Applying this to ()7.1|) it follows from f)7.5j) that the time averaged energy flux across the 
horizon is given by 

^(£ in - E out ) = l -A <$> 2 u(u - mQ h ) , (7.7) 

where A is the proper "area" of the horizon. Thus, in the frequency range of equation (|7.2jl . 
the energy loss across the horizon is negative and hence superradiance is obtained. 

8 Conclusions 

In this paper we have studied the phenomena related to the acoustic geometry of the axially 
symmetric accretion flow. We have discussed the formation of the shock and acoustic horizons 
using a general relativistic formalism in a model of the flow based on a perfect polytropic 
gas. We have found that a multitransonic flow can form with two sonic points and a shock 
in between. Null curves depicted in figures IU1 and 1101 clearly demonstrate the presence of the 
acoustic black hole at regular sonic points and of the white hole at the shock. 

Multitransonic accretion happens for a range of values of the specific energy density not 
much exceeding the rest mass. The analogue Hawking temperature of the inner horizon, 
calculated for a typical set of physical parameters, is lower but close in magnitude to the 
black-hole Hawking temperature. However, at the middle horizon that forms at the shock 
location, the surface gravity is formally infinite owing to discontinuities in the speed of sound 
and the radial velocity. We expect that in a realistic fluid with viscosity, discontinuity will be 
removed and replaced by a steep slope resulting in a finite but possibly very high analogue 
Hawking temperature. 
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